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ABSTRACT 

We have used the new extended A configuration of the IRAM Plateau de Bure interferometer to study the dense molecular gas in 
the nucleus of the nearby spiral galaxy NGC6946 at unprecedented spatial resolution in the HCN(l-O) and 12 CO(2-l) lines. The gas 
distribution in the central 50 pc has been resolved and is consistent with a gas ring or spiral driven by the inner 400 pc long stellar bar. 
For the first time, it is possible to directly compare the location of (dense) giant molecular clouds with that of (optically) visible HII 
regions in space-based images. We use the 3 mm continuum and the HCN emission to estimate in the central 50 pc the star formation 
rates in young clusters that are still embedded in their parent clouds and hence are missed in optical and near-IR surveys of star 
formation. The amount of embedded star formation is about 1.6 times as high as that measured from HII regions alone, and appears 
roughly evenly split between ongoing dust-obscured star formation and very young giant molecular cloud cores that are just beginning 
to form stars. The build-up of central mass seems to have continued over the past > 10 Myrs, to have occurred in an extended (albeit 
small) volume around the nucleus, and to be closely related to the presence of an inner bar. 
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1. Introduction 

The nearby (D=5.5Mpc, l"=27pc) Scd spiral galaxy NGC 6946 
is an ideal laboratory to study the causes (and consequences) 
of intense star formation in the vicinity of a galaxy nucleus, al- 
though its low galactic latitude hampers studies at optical and 
shorter wavelengths. Within its central 300 pc, NGC 6946 is cur- 
rently undergoing an intense burst o f star formation, as indi- 
cated by strong far-infrared emission dDevereux & Youndl993l) 
and the presence of n umerous hydrogen recombination lines in 
its nuclear spectrum (Engelbrac ht et alj I1996I) . The proximity 
of NGC 6946, combined with recent technological advances in 
mm-interferometers which now routinely reach sub-arcsec reso- 
lution, enables a detailed study of the interplay between infalling 
dense molecular gas, active star formation, and the energetic 
feedback from young massive stars on scales of individual giant 
molecular clouds (GMCs). We have recently undertaken such a 
study which has shown that NGC 6946 appears to be a "textbook 
case" for molecular gas resp onding to the gravitat ional potential 
of a small-scale stellar bar (ISchinnerer et al. 2006, hereafter pa- 
per I). The gas flows inward along an S-shaped spiral structure 
and accumulates in a massive nuclear clump with a size of about 
2" (60 pc). This clump contains about 1.6 x 10 7 M of molec- 
ular gas, and appears to completely obscure the very center of 
NGC 6946. 



* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). 



In this letter, we describe results from recent nun- 
observations obtained with the new, expanded baselines of the 
IRAM Plateau de Bure interferometer (PdBI) which yield a spa- 
tial resolution of about 0.35", and thus allow a direct comparison 
with the best available (space-based) optical and/or near-infrared 
maps. The new PdBI observations were designed to address a 
number of open issues related to the connection between active 
star formation and the gas properties. Of particular interest is the 
question how the sites of active star formation compare to the 
molecular gas flow, and whether optical or even near-infrared re- 
combination lines yield a complete picture of the current star for- 
mation rate in the nucleus of NGC 6946. We briefly describe the 
observations and data reduction procedures in §|2] and present 
the resulting maps of the 12 CO(2-l) and HCN emission in §[3] 
In §|4] we analyze the gas dynamics and various star formation 
tracers, and briefly discuss the implications of our results. 

2. Observations and Data Reduction 

We used the dual-receiver capability of the PdBI to obtain si- 
multaneous observations of the HCN(l-O) and I2 CO(2-l) lines 
at 3 mm and 1 mm, respectively. The 9 hr long observations were 
performed on January 15th, 2006, with all six PdBI antennas in 
their new, extended A configuration providing baselines between 
136 m and 760 m. For calibration and mapping, we used the stan- 
dard IRAM GILDAS softw are packages CLIC and MAPPING 
dGuilloteau & Lucai l2000h . The phase center of the observa- 
tions was at 20h34m52.36s +60d09ml5.96s (J2000.0). All ve- 
locities are reported relative to v sys (LSR) = 50 kms~'. During 
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Fig. 1. Intensity map of the 12 CO(2-l) line emission at 0.4"x0.3" 
resolution (color and contours). The black cross marks the posi- 
tion of the dynamical center as derived in paper I. The beam size 
and the linear scale are shown for reference in the bottom right 
corner. The contours start at 0.6 Jybearrr'kms -1 with the same 
step size. The red dashed ellipse indicates the location of the as- 
sumed bar (see paper I). 
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R.A. (J2000) 

Fig. 2. Comparison of the intensity map of the HCN(l-O) line 
emission (color) and the neighboring continuum emission (con- 
tours) at 1 .0"x0.7" resolution. The contours start at 3<t and are 
in steps of lcr=0.21 mJybeairT 1 . The black cross marks the po- 
sition of the dynamical center as derived in paper I. The beam 
size is shown in the bottom right corner. 



the observations, the typical system temperature in the sideband 
containing the line of interest at 89 (230) GHz was about 110 
(300) K. 

The new 12 CO(2-l) observations were combined with our 
previous dataset (presented in paper I) that was obtained with 
the same spectral setup in the old AB configuration. The data re- 
duction and imaging closely followed the procedures described 
in paper I. In order to obtain emission line-only maps, we 
subtracted (in the uv plane) continuum maps obtained from 
line-free channels over the following spectral ranges: 88.335- 
88.579GHz, 88.654-88.704GHz, and 88.767-88. 886GHz for 
the 3 mm band, and 230.21 1-230.367 GHz and 230.605- 
230.761GHz for the 1 mm band. Using uniform weighting, 
the resulting CLEAN beams are 0.97" x 0.65" (PA 93°) and 
0.40" x 0.29" (PA 90°) for the HCN(l-O) and 12 CO(2-l) line, 
respectively. For both lines, we defined a CLEAN region based 
on the 0th moment map. The rms in the 6 km s _1 wide channels 
is 2.3 mJy beam- 1 for HCN(l-O) and 4.8 mJy beam-' for CO(2- 
1). The rms in the continuum maps is 0.2 (0.21) mJybeanT 1 
and 0.7 mJybeairT 1 for natural (uniform) weighting at 3 mm 
and 1 mm, respectively. Moment maps were calculated using 
the GIPSY task 'MOMENTS' requiring that emission is present 
above a 3(2)<x clipping level in at least two adjacent channels of 
the CO (HCN) data cube. 



3. Results 

3.1. Distribution and Kinematics of CO 

As evident from Fig.Q] the new 12 CO(2-l) observations clearly 
resolve the structure of the inner CO spiraQ described in paper 
I. The south-eastern gas lane splits into two components while 
the north-western lane still appears as a single structure. More 
importantly, the structure of the nuclear clump is resolved into 
two peaks north and south of our derived dynamical center with 
extensions to the east and west which connect to the gas lanes. 
Overall, the molecular gas is in clumps that are mostly unre- 
solved at our resolution of (1 lx8)pc 2 , suggesting that the GMCs 
are rather compact, dense, and gravitationally bound. 

The emission lines in the central ~ 3" have complex profiles 
we fit with 2- or 3-component Gaussians. They often have two 
peaks, at (-45 + 10)kms-' and (+35 + 10)kms _1 relative to the 
systemic velocity, with a third peak at (-5 +5)kms _1 in the inner 
0.4", near the dynamical center (Fig. [3]). The velocity dispersions 
are 20 to 30 km/s in the two main peaks, and < 10 km s _1 in the 
third, central component. Most of the negative velocity compo- 
nent is confined to one side of the dynamical minor axis, a result 
which is consistent with the gas kinematics being dominated by 
circular motions. This is also true for the positive velocity com- 
ponent besides an extra peak about 0.3" NE of the dynamical 
center, close to the northern CO peak in the integrated map. This 
is also illustrated by the double-peaked line profile of the north- 
ern peak (Fig. [3]) when compared with the simpler line profile of 
the southern peak. 



1 For the description of the various morphological components, we 
adopt the nomenclature introduced in paper I. 
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Fig. 3. Gas kinematics in the central 3" of NGC 6946. Middle 
panel: Integrated flux of the negative (positive) velocity com- 
ponent from a double Gaussian fit [see text; white (black) con- 
tours] overlaid on the total integrated 12 CO(2-l) flux (color). 
Top/Bottom panels: Typical 12 CO(2-l) spectra (from single pix- 
els) and their corresponding Gaussian fits display a double- 
peaked profile throughout the nuclear region; the profile close 
to the dynamic center shows a third component close to the sys- 
temic velocity. The velocity is given relative to the systemic ve- 
locity and the y-axis is line intensity in JybeairT 1 (with a tem- 
perature conversion of 200 KJy -1 ). 



3.2. HCN(1-0) Emission 

The HCN(l-O) emission associated with the nuclear and the 
southern clump is shown in Fig. [2] The map contains a total flux 
of about 4.5 Jy kms~' above a 2cr threshold, that is about 8% 
of the flux of 11.4K kms " 1 measured for the central position 
bv lGao & Solomon! d2004l) using the IRAM 30m telescope (as- 
suming a temperature-to-flux conversion of 4.8 Jy Kr 1 ). 60% of 
the recovered flux originates in the nuclear clump. The average 
fine intensity ratio of HCN(1-0)/ 12 CO(2-1) is 0.02-0.03. Thus, 
the non-detection of the northwestern part of the spiral structure 
is probably due to our S/N limit combined with the lack of short 
spacings in the HCN data. While there is a good correspondence 
between the HCN and 12 CO(2-l) emission in the southeastern 
spiral structure, the integrated HCN emission has an additional 
peak in the nuclear disk about ~ 1" west of the dynamical center. 



3.3. Millimeter Continuum 

Continuum emission has been detected at 3 mm and, tentatively, 
also at 1 mm. The integrated 3 mm flux is ~ 2.8 mJy, with a peak 
flux of 1.9 (1.7)mJybeam _1 corresponding to a 9(8)<x detection 
with natural (uniform) weighting. The position of the peak is at 
20:34:52.307 +60:09:14.37 (J2000) and offset by ~ 0.4" from 
the dynamical center as derived in paper I. At the low levels 
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Fig. 4. Top: HCN(l-O) intensity map (contours) overlaid on the 
12 CO(2-l) one (color). The HCN contours start at 50% of the 
maximum in steps of 10%. Both beams are shown in the bottom 
right corner with the HCN one being the larger one. Bottom: 
HST Paa line emission (contours) overlaid on the VLA 6cm 
radio continuum (color). The position of the dynamical center 
(cross) and the 3mm continuum peak (star) are indicated. The 
6cm beam is shown in the bottom right corner. 



the 3 mm continuum shows an extension to the west (Fig. [2]) 
and is slightly resolved in the uniformly weighted map with a 
Gauss-fit of size (FWHM) 1.22" x 1.06" (PA ~ 96°). Continuum 
emission at 1 mm is present at the 4cr level with a peak flux 
of 3.0 mJy beam -1 in the naturally weighted data. A recent 0.3" 
resolution 6 cm continuum map bv lTsai et alj d2006) shows that 
the nuclear radio continuum breaks up into two bright peaks with 
extensions to the north and west. The 3 mm continuum peak is 
off-set by ~ 0.4" to the north of the bright 6 cm peak (Fig. |4| 
strongly indicating that the 3 mm and 6 cm continuum emission 
have a different origin. 



4. The Central 100 pc: Gas Flow and Star Formation 

4.1. Bar-Driven Gas Kinematics 

In paper I, we showed that the CO spiral can be explained by the 
response of the gas flow to the gravitational potential of the in- 
ner NIR bar in the central 300 pc. Our qualitative model requires 
the existence of an Inner Lindblad resonance (ILR) at a radius 
of ~ 2.2" due to a -400 pc long secondary bar. In this context, 
we expect gas to accumulate at or close to the ILR and, thus, star 
formation to occur. Such a mechanism is evoked to explain the 
formation of the nuclear star formation rings of kiloparsec di- 
amete r in barred galaxies (e.g. IWada & Habell99"2l:lAllard et alj 
2006). While the CO data of paper I showed streaming motions 



4 



Schinnerer et al.: Central 50 pc of NGC 6946 



associated with the bar, the exact distribution of the gas in the 
nuclear clump could not be resolved. 

The CO(2-l) intensity distribution is very reminiscent of the 
so-called 'twin peaks' struct ure seen in lower re solution CO(l- 

0) images of barred galaxies dKennev et aljl9 92) and associated 
with gas accumulation at the ILR. The two peaks of the nega- 
tive and positive velocity component are located at a radius of 
~ 0.5" (Fig. |3), i.e. within the ILR of our model. Dense gas as 
traced by the HCN(l-O) line roughly coincides with those CO(2- 

1) peaks (Fig. |4j». One way to interpret these gas peaks is as the 
contact points of two spiral arms connecting in the central lOOpc. 
Our derived line intensity ratios (using the integrated flux and 
correcting for different beam sizes) correspond to temperature 
ratios of fjggffzjU ~ 0.15 - 0.20, similar to values found in kpc- 
sized starburst rings (e.g. lKohnoet~aT1ll999l) . As the presence 
and distribution of dense gas is analogous to what is observed 
in kpc-sized starburst rings of large-scale bars, our new data to- 
gether with the kinematic analysis (presented in detail in paper 
I) are consistent with the interpretation that the inner bar drives 
the CO spiral even within the nuclear clump, 

The nature of the additional eastern peak in the positive ve- 
locity component of the CO(2-l) emission is not clear. However, 
the young massive stars present there (Fig. @|i will clearly alter 
the properties of the (dense) molecular gas. In addition, we start 
to probe size scales where the size and stability of the GMCs 
themselves is becoming important. Observations of other molec- 
ular gas tracers might therefore be needed to clarify the situation. 



about 0.1 (within the uncertainties on the 1mm flux), which is 
expected f or free-free em ission. Following the equations sum- 
marized in lJohnsorJ d2004l) we find an ionizing luminosity Qlyc 
of ~ 10 52 s~', equivalent to a SFR of 0.1 M yr~'. This is an 
additional contribution to the total SFR, originating in regions 
which are n ot (yet) seen in hydroge n recombination lines. 

Finally. iGao & Solo mon (2004) showed that there is a lin- 
ear relation between the luminosity of the HCN(l-O) line and 
the far-infrared luminosity. This relation appears to not only 
hold true for entire galax ies, but also fo r individual star forming 
cores within our Galaxy dWu et al.l2005l) . Thus, it seems reason- 
able to use the HCN(l-O) line flux as a proxy to estimate the 
rate of stars that are in the process of forming. Our measured 
HCN(l-O) line flux of 2.5 Jykms 1 corresponds to a luminos- 
ity of L'(HCN)=3.1 x 10 5 Kkms-'pc 2 . Following eq. [11] of 
IGao & Solomon] (|2004), we obtain another contribution of 0.06 
M Q yr" 1 to the total SFR. 

In summary, the different diagnostics (Paa/6cm, 3mm, HCN) 
all trace star formation, but in different evolutionary stages. This 
is also apparent from their different spatial distributions. We 
find about twice as much embedded star formation (traced by 
HCN and 3 mm continuum) as star formation which has already 
emerged from its dust cocoon (traced by Paa). This implies that 
there has been significant massive star formation in the central 
lOOpc of NGC 6946 over the past >10 Myr. The 'older' star for- 
mation has neither disrupted the GMCs seen in HCN nor has it 
prevented the nuclear star formation evident in the 3 mm contin- 
uum. 



4.2. Census of Nuclear Star Formation 

The high extinction present in the central lOOpc makes it diffi- 
cult to estimate the true star formation rate (SFR). In what fol- 
lows, we compare different SFR indicators in the central 3"x3" 
to develop a sense for how much star formation is still deeply 
embedded in dense gas clouds and hence does not reveal itself 
in observable emission from hydrogen recombination lines. 

The 'visible' SFR can be traced by the NIR Paa line emis- 
sion arising from HII regions. Us ing a Paa flux of 2.9 x 
10~ 13 erg s _1 crrT 2 from the image by Boker et al. ( 1999) a nd as - 
suming an extinction of Ay — 4.6mag (Ouillen & Yukita 2001), 
we derive a SFR of 0.13 M yr~' (see paper I for details). 
This value can be directly compared to the SFR rate derived 
from the non-thermal radio continuum which is associ ated with 
synch r otron emission from supernovae explosions (e.g. lCondonl 
119921) . iTurner & Hoi dl983l) find that most (~ 90%) emission at 
6cm in the nuclear region of NGC 6946 is non -thermal in origin, 
while the detailed work by Tsai et alj d2006l) show that point- 
like sources (either due to HII regions or supernova remnants) 
contribute about 30% to the total 6cm flux. We measure a flux 
densit y of abou t 15 m Jy in the central lOOpc. Using equation 
T211 of ICondonld 19921) and assuming that 80% is non-thermal in 
origin, we derive a SFR of 0. 1 M Q yr~' (also taking into account 
stars below 5 M ) in agreement with the SFR derived from the 
Paa line. This agreement might be expected, as the spatial distri- 
butions of Paa and 6 cm radio continuum agree very well (Fig. 
@), and it supports our assumption that both are tracers of stellar 
populations of similar age. 

The millimeter continuum emission, on the other hand, 
traces free-free emission from star formation sites that are still 
embedded in their parent clouds. This picture is supported by 
two arguments: first, the distribution of the mm continuum emis- 
sion differs from that of the radio continuum. Secondly, the spec- 
tral index aj,-i between the 3 mm and 1 mm continuum flux is 



4.3. Nuclear Cluster Formation? 

The nucleus of most late-ty pe spiral galaxies is marked by a 
compact nuclear star cluster dBoker et alj|2002l) . In NGC 6946, 
however, such a nuclear cluster cannot be identified. Possibly, 
a nuclear cluster exists, but is obscured by the large amounts 
of molecular gas and dust in the central 1". Another possibil- 
ity is that we currently witness the birth of a nuclear cluster. 
According to our model (paper I), the bar-driven CO spiral struc- 
ture requires the presence of an ILR, which itself implies a mass 
concentration on scales of about 20-30 pc (~ 1"). Comparison 
of the different star formation tracers shows that stars are being 
formed in the central 60 pc at a rate of about 0. 1 M yr~ 1 over the 
past > 10 Myrs translating into a total of roughly 10 6 M in stel- 
lar mass. The spatial distribution of the SFR tracers shows that 
star formation occurs in distinct areas (clusters, GMCs) spread 
over the central 2" (60 pc). This suggests that the build-up of 
central mass takes place in an extended (albeit small) volume 
around the nucleus, rather than at the very center. This is consis- 
tent with the recent identificati on of blue, extend ed disks around 
a number of nuclear clusters dSeth et al.ll2006l) . NGC 6946 of- 
fers the unique opportunity to study the processes linked to the 
build-up of a central mass concentration. 
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